Fluctuation-Driven First-Order Transition in Pauli-limited <i-wave Superconductors 
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We study the phase transition between the normal and non-uniform (Fulde-Ferrell-Larkin- 
Ovchinnikov) superconducting state in quasi two-dimensional d-wave superconductors at finite tem- 
perature. We obtain an appropriate Ginzburg-Landau theory for this transition, in which the fluctu- 
ation spectrum of the order parameter has a set of minima at non-zero momenta. The momentum 
shell renormalization group procedure combined with e expansion is then applied to analyze the phase 
structure of the theory. We find that all fixed points have more than one relevant directions, indicating 
the transition is of the fluctuation-driven first order type for this universality class. 
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It was pointed out forty years ago [1,2] that a super- 
conducting state with an inhomogeneous order parameter 
can be stabilized by a Zeeman splitting between electrons 
with opposite spins, that is comparable to the energy gap. 
This inhomogeneous superconducting state, or Fuldc- 
Ferrell-Larkin-Ovchinnikov (FFLO) state, has been the 
subject of mostly theoretical study for many years [3]. 
The situation has changed recently as experimental re- 
sults suggestive of the FFLO state start to emerge [4-11]. 
We would like to mention that recent experimental re- 
sults on the heavy fermion compound CeCoIns, a quasi 
two-dimensional (2D) d-wave superconductor, are partic- 
ularly compelling [9-11]. It is worth pointing out that the 
FFLO state may also be realized in high density quark 
matter, and is thus of interest to the particle physics 
community [3]. 

Given the long history of the subject matter, it is per- 
haps somewhat surprising that thus far most of the the- 
oretical studies of the FFLO state are of the mean-field 
type. On the other hand one expects that quantum and 
thermal fluctuation effects should be much more signif- 
icant in FFLO superconductors than in ordinary BCS 
superconductors, as the FFLO phase breaks the transla- 
tional symmetry in addition to the gauge symmetry. In 
a recent work [12], one of us studied the FFLO phase 
in quasi ID superconductors, and used bosonization to 
treat intra-chain quantum fluctuations exactly; one of 
the conclusions was that the transition from the FFLO 
phase to the BCS phase is a continuous transition of the 
commensurate-incommensurate type, in contrast to the 
first order transition commonly asserted in the literature. 
The effect of thermal fluctuations of the superconduct- 
ing order parameter was discussed by Shimahara [13]; he 
argues that the enhanced fluctuation effects destabilize 
certain types of mean-field FFLO states. Anisotropy in 
pairing or electron dispersion may suppress these fluctu- 
ation effects however [13]. 

It is perhaps not quite well recognized yet that the 
FFLO superconductors are realizations of the Brazovskii 
model [14], which describes a large class of statistical me- 
chanical systems in which the fluctuation spectra of the 



order parameter have their minima away from the ori- 
gin in momentum space; this is precisely the case for the 
FFLO state, which prefers the superconducting order pa- 
rameter to carry finite momenta. In its original form, the 
Brazovskii model assumes that the fluctuation spectrum 
has a continuous set of degenerate minima; it was shown 
that the fluctuation effects are so strong that they render 
the transition between the ordered and disordered phases 
a fluctuation-driven first order transition, even if the 
mean-field theory suggests a second-order transition. In- 
terestingly, the transition between the normal and (possi- 
bly) FFLO phase in CeCoIns was indeed found to be first- 
order [9,10]. There are two possible origins for the first 
order nature of the transition. Firstly, it is known that 
near the tricritical point where the normal, FFLO and 
BCS superconducting phases merge within mean field 
theory (at T as 0.56T C ), the effective Ginzburg-Landau 
free energy has a contribution from quartic terms that 
is in certain cases negative [15,16], in which case the 
mean-field theory itself predicts a first-order transition. 
In CeCoIns however, the FFLO phase was observed only 
at temperatures much lower than the tricritical point. 
It has been pointed out recently [17] that in the low- 
temperature regime the quartic term makes positive con- 
tributions to the free energy; thus a second-order transi- 
tion would be expected at mean-field level. If this is the 
case then the origin of the first-order transition must be 
due to fluctuation effects left out at the mean-field level. 
It is this second possibility that we focus on in this work. 

While the FFLO phase and the original Brazovskii 
model share the common feature that the order param- 
eter spectrum has it minima away from the origin, one 
cannot apply the Brazovskii results directly to real sys- 
tems like CeCoIns because in such systems there is al- 
ways anisotropy in either pairing potential or electron 
dispersion, which would generically reduce the continu- 
ous set of degenerate minima to a discrete set. Thus in 
this work we take this important effect into account, and 
use a model that contains the anisotropy that is appropri- 
ate for a d-wave superconductor with four-fold symme- 
try. We perform a renormalization group (RG) analysis 
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(combined with an appropriate e expansion, see below), 
and show that the transition is generically first order (at 
least when e is sufficiently small), even when the mean- 
field analysis suggests a second-order transition. In the 
following we first outline the derivation of an appropri- 
ate Ginzburg-Landau theory from a microscopic pairing 
model, then perform the RG analysis. 

We consider a weak-coupling quasi 2D d-wave 
superconductor, whose partition function is Z = 
J D[&\ *] exp{-S}, with S = S + S int , where 

So=J2 [«^n-C(k)-<^*£(k,w n )*a(k,w„), 
Sint = -T E ^,k'*j(k + * v + w ) 

k,k',q ,u> 

■*{(-k, -i/)* T (-k', -i/)*i(k' + q, !/' + w), (1) 

with \I>,\I>t being Grassman variables. Index cr =f, J, enu- 
merates electron spin, and I is the Zeeman splitting that 
stabilizes the FFLO state [18]. The electron dispersion 
£(k) = e(ky) + J cos k z d — ep, where ky denotes the mo- 
mentum parallel to the planes, J is the strength of hop- 
ping between the layers. Interaction Vk.k' = V /k/k with 
V > 0, is assumed for simplicity to be independent on 
the z-component of momenta with /k = cos2#k for the 
d-wave pairing. We consider here the clean system only, 
assuming that disorder does not affect qualitatively the 
phenomena under consideration. 

We decouple then the quartic term via the Hubbard- 
Stratanovich transformation by introducing the complex 
field A(q, u n ), that serves as a superconducting order 
parameter. To obtain the Ginzburg-Landau (GL) ac- 
tion in powers of A, we perform subsequently the cumu- 
lant expansion integrating out the fermions having the 
Green function (<r = ± for up and down spins respec- 
tively) G _1 (q, u n ) = iui n — £(q) — al. Considering here 
the finite-T transition only, we retain the zero-frequency 
component A(q, oj n = 0) in all cumulants. The resulting 
functional takes the form 

^ = 5> 2 (q)|A(q)| 2 
q 

+ V X4(qi,..,q4)A( qi )A(q 2 )A*(q3)A*(q4), (2) 

* — 'qi,..,q 4 

where the prime in the sum over momenta in the quartic 
term indicates that the condition qi + q 2 — q 3 — q 4 = is 
taken into account. Since the transition occurs into the 
state that is non-uniform in space, the full momentum de- 
pendence of K 2 and K4 must be kept [19]. K 2 (q) is given 
by the standard bubble diagram with two external legs, 
K 2 (q) = l/V-T£ ki „ B /kG| (k + q,i/ n )xG ; (-k, -u n ). 
In this formula, one performs then straightforwardly the 
summation over frequency and integration over momen- 
tum within the shell around the Fermi surface \k~kp\ < 
wd/vfQ<-f)- The form of the Fermi surface is assumed 



to have the same 4-fold d-wave symmetry in the fcii 
plane. Distinguishing also the components of the Fermi- 
momentum parallel and perpendicular to the planes, 
hp = (k|| F \ k z F ^ ) , Vf = (v||,w z ), we find as a result 
that for q < hp [20] 

1 1 /■*■/<* k\pdk z r 

In the equation above, h = 1, 

z q = \ \ v \\% cos (#v - 0q) - Jdq z - sin k z d] + I, (4) 

with 8, 8 V and q being the in-plane angles of k, vp and 
the pairing momentum q respectively, k^, vp as well as 
8 V , are themselves functions of k z and 8 characterizing 
the Fermi surface. 

To determine the absolute value of the pairing momen- 
tum g and its direction, it is necessary to find the minima 
of K 2 (q) with respect to gy and # q , as well as q z . It is 
clear that the ordering wave vector must lie in the (q x , q y ) 
plane meaning that q 0z = 0. However, investigation of 
_ftT 2 (q) regarding the minimum with respect to # q , re- 
veals that Eq. (3) has extrema for the two sets of values: 

Oi = f,x>X>-f and ^ = O'f.^TT correspond to 
the nodal and anti-nodal directions in the {q x , q y ) plane. 
Generally speaking, for each of these sets one obtains 
the different values of qo as a result of the solution of 
equation dK 2 (q) /dqo . Whether both of these sets mini- 
mize K 2 (q), or only one of them is actually a minimum 
with the other being the maximum, is determined by the 
specific form of the Fermi surface. It is important that, 
if both sets with the corresponding values of q are the 
minima, the actual transition will occur into configura- 
tion described by the set leading to the largest critical 
value T c [20,21]. In case of the anti-nodal ordering, it 
is possible to expand K 2 (q) for q located in the pockets 
near the minima ?o2)<?g ' corres P on ding to the direction 
8 U K 2 {q) = r + ai l \q x - g<2) 2 + a« \q y gg) 2 + 79 2 , 

where a x ^ and 0$ mutually interchange for the neigh- 
boring pockets. For the case of nodal ordering, one can 
show from Eq. (3) that K 2 (q) = r + a(q x - g^) 2 + 

2/3 (2) fe - q£l){q y - <?g) + a(q y - gg) 2 +7?* , where /?« 
are opposite in sign for the neighboring pockets. By the 
simple rotation of coordinate system in the (q x ,q y ) plane 
by 7r/4, the latter expansion, however, reduces to that 
for the anti-nodal case. 

The quartic kernel in Eq. (2) is given by the bubble 
containing four electron Green functions and four exter- 
nal legs representing the order parameter field A(q), 

K 4 = (T/4) £ {/^ +q4 _ qa Gt(k + qj ) i/ n )G i (-k,-i/ n ) 

^n,k 
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G|(k + q 4 , j/„)G|(— k + q 2 - qi,-v n ) + [q4 -> q3]} , (5) 

where [q 4 — > q 3 ] stands for the same expression as right 
before, with only q 4 substituted by q 3 . The kinematic 
constraint qi +q 2 = q3 + q 4 is implied in Eq. (5). It will 
not be required, however, to know this cumulant for all 
values of momentum variables. Since we are considering 
the renormalization group treatment involving the wave- 
vectors located in the pockets near q^ , only those values 
in Eq. (5) are of interest, in which q's point right to 
the centers of the pockets and satisfy the aforementioned 
constraint . The following distinct possibilities can be 
readily enumerated, once one denotes by i the number of 
the pocket in the qy plane; i equals to 1,2,3,4 starting 
from that with the lowest value of angle, with the formal 
condition i + 4 = i. 



^(q^.q^.q^.qS ) = «o/4, i = 1, ...,4; 



(6) 



^4(q ,Qo iQo >Qo J 



V2 /4,i=l,..,4; (7) 



7v- (_(*) _(*+2) _(*) _(*+2)\ _ „ M ,• _ 1 9. 



^4(q ,Qo )Qo '^o J 



u/4. 



(8) 



(9) 



Looking at the expressions given by Eqs. (6)- (9), we 
see that uo describes the interactions between the modes 
within the same pocket, while the other parameters ac- 
count for the inter-pocket scattering. In parts with 
and u„, the interaction occurs between the pockets that 
have the angles between their qo's equal to 7r/2 and it 
respectively. Without providing the explicit expressions 
for those coefficients, we note that the only point impor- 
tant in the general derivation is that all interactions are 
non-singular at the critical values T c and qo- The in- 
teractions can in principle have arbitrary signs that may 
change along the critical line T c = T C (I) on the (T, I) 
plane. For example, as it was mentioned in the intro- 
duction, interaction w is negative close to the tricritical 
point, where the normal, uniform and non-uniform super- 
conducting phases meet [16]. At the same time, at lower 
temperatures uq seems to be positive [3]. If u n < 0, the 
transition is necessarily first order already at the mean- 
field level. Hence, we will assume that we consider the 
transition only in those regions on the phase diagram, in 
which at least uq is positive. 

To distinguish the modes with the wave vectors belong- 
ing to the different pockets, it is convenient to introduce 
the shifted momenta k = q — q 1 ' and decompose the to- 
tal field A(q) into the parts A,(k = q — qg ). Each part 
Aj(k) accounts for the fluctuations having the momenta 
in the vicinity of q ^ . It is clear that under such decom- 
position, the kinematic constraint ki + k 2 = k 3 + k 4 for 
the shifted vectors is preserved. Though kj in the arising 



quartic terms generally belong to the different pockets, 
one can treat them during the formulation of RG equa- 
tions, as if they are located in one and the same pocket 
around the origin. We will use in RG equations below the 
form of the propagator obtained for the anti-nodal or- 
dering assuming for clarity that = a.\ ^ = a 2 , 
meaning the spatial anisotropy in the spectrum of ex- 
citations. The issue of spatial anisotropy in RG near 
quantum critical points was addressed in different phys- 
ical context in Ref. [22], albeit the anisotropy there was 
related to fermionic excitations. 

After the appropriate rescaling of momentum variables 
and fields, the general GL action takes the form: 

T = E E [ r + a * )k * + a v )k l + k *} i A *( k )i 2 

i k 

+ E' {fe;} | (W 4 ) E A i (k 1 )A i (k 2 )A*(k 3 )A*(k 4 ) 
+^ /2 ^A,(k 1 )A, +1 (k 2 )A*(k 3 )A* +1 (k 4 ) 

+«, J2 A i (k 1 )A i+2 (k 2 )A*(k 3 )A* +2 (k 4 ) 

+v (A 1 (k 1 )A 3 (k 2 )A*(k 3 )A:(k 4 ) + c.c)} . (10) 

A few more remarks on the notations in Eq. (10) are in 

order. The notation } implies that the summation 
over kj is taken with the restriction ki + k 2 = k 3 + k 4 . 

means that the sum over i is performed with the 
condition i + 4 = i. In all terms of the quartic part, 
except that with wo/4, the permutational symmetry be- 
tween the fields arising from the obvious permutations of 
momenta in arguments of Eqs. (7)-(9), is taken care of 
explicitly, canceling thus the factor of 4 in denominator. 
This greatly simplifies the subsequent RG loop analysis. 

Couplings uo and « T are in fact the primary parame- 
ters, whose flow under rescaling determines the character 
of transition. To see this, we calculate the free energy 
at the mean-field level for two possible phases: l)Fulde- 
Ferell (FF) phase with A(r) = A e lqr and 2)Larkin- 
Ovchinnikov phase having A(r) = Ao cos(q • r). The 
values are J^ff = —\r\ 2 /u , .Flo = -2|r| 2 /(u + 2u v ). 
The LO phase has the lower energy when u > 2u v , 
while the FF phase is more favorable under the opposite 
condition. The considerations above necessarily imply 
that u > 0, since only in this case the transition is of 
the second order at the mean field level. In addition, 
if LO phase is realized, one must ensure that not only 
u > but also u + 2u 7T > 0. Those requirements will 
be presumed fulfilled in the subsequent treatment. 

Simple tree-level scaling applied to Eq. (10) shows that 
if the effective dimensionality of the problem, d > d c = 4, 
the interactions are irrelevant and the transition is of 
the second order. To proceed, we will generalize the 
z-component of momentum to k± having the dimen- 
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sionality 2 — e, and integrate out the modes in the 
thin layer A/e l < k x ,k y < A around the square shell: 
-A < k x < A, k y = ±A; -A < k y < A, k x = ±A, with 
the integrals over k± taken from — oo to oo. The arising 
one-loop RG equations for interactions are: 

duo/dl = eu - 2/ [(5/4)u§ + 2ul /2 + u*] , (11) 



du^/dl = eu^ -2f uqUk + u 2 + u 2 /2 + v 2 /2 , (12) 



du 



n/2 



dl 



with 



= - 2 



fu n / 2 (u + u n ) + g(ul /2 + — ) 



(13) 



dv/dl = ev - 2v [fu n + 2gu n / 2 ] 



(14) 



(2tt)3 



k±dk± 



' ,A 

/ dk x g 2 { ai kl,a 2 K 2 ) 
Jo 



f 

Jo 



dk y g' z (a 1 K 2 ,a 2 k 2 v ) 



the fixed points in the whole range < a < 1. We thus 
find no stable fixed points at the one-loop level, meaning 
that one needs to tune at least two parameters (r and v) 
to reach the fixed points. This implies that the transition 
will be generically of the first order, even if the mean-field 
theory suggests a second-order transition. The situation 
here is similar to that near transitions described by the 
effective Hamiltonians of anisotropic systems [23]. 

In summary, we have obtained an effective Ginsburg- 
Landau theory for the transition from normal to the 
FFLO state in quasi 2D d-wave superconductors at all 
non-zero temperatures. RG analysis of the theory indi- 
cates that the transition is generically first order, even 
when the mean-field theory suggest a continuous transi- 
tion. This fluctuation-driven first order transition is due 
to the enhanced fluctuations of the FFLO state, associ- 
ated with additional broken symmetries. Our result is 
consistent with the first order character of the transition 
observed in CeCoIns. 
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and the NSF Grants No. DMR-0225698 (KY) and PHY- 
9907949 (at KITP). 
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/oo [>A 
k±dk± / dk x Q(a\k 2 ,a 2 K 2 
oo Jo 



,A 

Q{a 2 k 2 x ,aik 2 ) + / dk y g(a 1 A 2 ,a 2 k 2 )g(a 2 A 2 ,a 1 k 2 ) 
Jo 



In the equations above, Q(x, y) = l/(kj_ + x + y), and we 
set r(l) <~ 0(e) to zero in Q(x, y) in the one-loop approx- 
imation. Looking for the fixed points, we absorb / by 
rescaling the interactions, generating thus the anisotropy 
parameter a = gj f = (2y/rj/(rj—l)) arcsin(?7— l)/(r] + 1), 
where rj = TaeiK{ai/a 2 , a 2 /ai} > 1. This parameter is 
cutoff- independent and satisfies < a < 1. As it can 
be seen from Eq. (14), it is reasonable to search sep- 
arately the fixed points that have v* = and v* ^ 0. 
We have found that setting the condition v* ^ 0, leads 
to the absence of fixed points in the space of real vari- 
ables for < a < 1. Concerning the fixed points having 
v* = 0, it can be shown that apart from the completely 
unstable Gaussian fixed point one has four more points: 



U n/2 = U 



= u* = 0,«S = 2e/5; 



V/2 



= 0,< = u* /2 = e/6; 
plus two more fixed points that are some cumbersome 
functions of a not to be presented here. For a = 1, one 
easily finds the latter to be M*y 2 = w* = Uq/2 = e/8, 
u t/2 = u n ~ M o/6 = e/16; while for a — > they both 
collapse onto the point m*, 2 = 0, u* = u* Q /2 — e/6. All 
of the found this way fixed points are unstable, since, as 
follows from Eq. (14), there will be at least one direction 



with the eigenvalue A„ = e — Au 



Tt/2 



2u*, positive at all 
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